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Abstract:  The  need  exists  for  robust  and  efficient  optimal  design  methods  for  application  to 
multibody  systems,  in  which  the  components  to  be  designed  represent  connections  between  large 
displacement,  large  rotation  motions  of  the  subsystems’  bodies.  A  specific  application  is  an 
occupant  restraint  systems,  such  as  the  Gunner  Restraint  System  (GRS),  in  which  both  the  vehicle 
and  the  gunner  can  undergo  large  relative  and  absolute  motions  under  extreme  driving  or  external 
threat  conditions.  In  addition,  the  restraint/connection  components  can  have  amplitude-dependent, 
time-dependent,  and  timing-dependent  behavior,  such  as  an  active  belt  retractor.  Current 
optimization  methodologies  are  ill-suited  for  this  problem,  suffering  from  infeasibility,  lack  of 
robustness,  and/or  high  computationally  expense.  This  paper  presents  an  extension  of  topology 
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optimization  techniques  to  consider  multibody  dynamics  systems  and  to  treat  the  much  more 
open  design  space,  which  can  include  passive,  active,  and  reactive  structures/devices.  The 
objective  is  to  obtain  an  optimally  combined  structural  and  material  system,  considering  the  best 
use  of  passive,  active  and  reactive  members.  This  paper  highlights:  1)  dealing  with  design 
objectives  that  consider  time -dependent,  dynamic,  large  deformation  responses;  2)  general 
representative  models  for  the  multi-disciplinary  (passive,  active  or  reactive)  components  in  a 
multibody  dynamics  simulation  system;  3)  designing  an  optimal  system  that  can  satisfy  multiple 
requirements  under  various  operating  conditions;  4)  an  efficient  sensitivity  analysis  method  for 
the  optimization  problem  of  the  restraint  system;  and  5)  a  general  and  advanced  optimization 
algorithm  that  can  solve  the  problems. 

Keywords:  topology  optimization,  multibody  dynamics,  sensitivity  analysis,  restraint  system, 
vehicle  safety,  automotive  vehicles,  active  devices. 

1.  Introduction 

Motivating  this  research  is  the  need  to  design  vehicle  occupant  restraint  systems  for 
improved  occupants’  safety  under  various  operating  conditions  and  often  hazardous 
environments.  Using  a  Gunner  Restraint  System  as  an  example,  the  occupant  (gunner) 
sits  or  stands  in  the  passenger  compartment  with  their  upper  torso,  arms,  and  head 
exposed  outside  the  top  of  the  vehicle.  The  restraint  system  should  not  only  be  able  to 
prevent  the  occupant  from  being  ejected  from  the  vehicle  but  also  be  able  to  assist  rapid 
entry  into  the  vehicle  during  a  rollover  or  other  accidents  to  avoid  injury  or  fatality.  For 
this  application,  the  restraint  system  should  also  help  stabilize  the  gunner  over  rough 
terrain  and  in  high  speed  maneuver  conditions  for  them  to  complete  their  functional  tasks. 
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The  restraint  system  may  involve  a  wide  range  of  possible  usage  of  passive,  active  and 
reactive  devices  which  could  be  mounted  at  many  possible  physical  locations  (interacting 
points)  between  the  vehicle  and  the  occupant.  These  devices  may  include  safety  elements 
such  as  belts,  airbags  and  retractors  and  may  have  to  be  activated  in  a  specific  sequence 
or  timing  to  protect  the  occupant  in  the  designed  situations.  For  the  purposes  of  this  paper, 
a  passive  device  is  defined  as  a  structure  or  device  that  responds  to  the  excitation 
passively  without  an  active  action.  An  active  device  is  defined  as  a  structure  or  device 
that  can  actively  respond  to  the  excitation  with  an  energy  supply  for  the  operation.  A 
reactive  structure  is  defined  as  a  class  of  smart  structure  that  can  react  to  external 
excitations  in  a  specially  designed  way  using  the  energy  pre-stored  in  the  system  or  from 
the  external  excitation  to  counteract  the  hazardous  loading  or  perform  other  desired  tasks. 
(Chiyo  et  al.,  2010,  Dong  et  al.,  2009;  Ma  et  al.,  2006a;  2007;  2008;  2010)  The  design  of  a 
restraint  system  must  also  consider  minimizing  the  system  weight,  complexity,  and  cost,  while 
maximizing  reliability,  durability,  and  occupant  friendly-ability. 

More  generally,  the  design  problem  of  interest  involves  multiple  multibody  dynamics 
systems  and  their  interconnections,  which  need  to  be  designed  to  constrain  the  relative 
motions/positions  of  the  multibody  dynamics  systems  for  given  objectives,  such  as  those 
related  to  the  safety  issues.  The  multibody  dynamics  systems  can  include  flexible  bodies; 
however,  in  this  paper,  we  limit  developments  to  rigid  multibody  dynamics  systems  for 
the  purpose  of  exposition.  The  application  focus  is  on  the  safety  system  design  problems 
related  to  automotive  vehicles,  including  military  vehicles,  such  as  gunner  restraint 
systems,  blast-protective  seating  systems  and  other  restraint  systems,  and  commercial 
applications,  such  as  passenger  safety  and  protection  systems  in  passenger  cars  for 
protection  against  crash  or  rollover.  Other  applications  vehicle  transportation  systems, 
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space  vehicle  landing  systems,  ground  and  sea  vehicles  mooring  systems.  For  a 
transportation  system,  the  design  objective  can  be  the  relative  movement  of  the  vehicle 
with  respect  to  the  carrier  vehicle  (ground,  sea  or  air)  for  a  transportation  task  in  a 
dynamic  environment.  The  design  space  could  include  connecting  chains,  networked 
belts,  or  other  constraint  mechanisms.  For  the  optimal  mooring  system,  the  design 
objective  could  be  the  vessel’s  lateral  and  longitudinal  accelerations  and  yawing 
movements.  The  design  space  can  be  all  the  possible  interactions  between  the  vessel  and 
the  dock  with  the  objective  to  find  the  optimal  mooring  system. 

Practical  solution  of  these  design  problems  requires  a  robust  and  efficient  optimal  design 
method  to  quickly  layout  an  optimal  restraint  system  between  the  multiple  multibody 
systems,  in  which  the  components  to  be  designed  can  represent  connections  between 
large  displacements,  large  rotation  motions  of  the  subsystems’  bodies.  In  addition,  the 
connection  components  can  have  amplitude-dependent,  time-dependent,  and  timing- 
dependent  behaviors,  such  as  that  with  an  active  belt  retractor.  Current  optimization 
methodologies  are  ill-suited  for  this  problem,  suffering  from  feasibility,  robustness, 
and/or  efficiency.  A  fundamental  multidisciplinary  structure  design  methodology  for 
multibody  dynamics  systems  is  presented.  This  design  methodology  identifies  optimally 
combined  multidisciplinary  structural  components  with  specific  geometric  and 
connectivity  configurations  and  also  mechanical  properties  for  the  given  (multiple) 
design  objectives.  One  challenge  in  developing  such  a  design  methodology  comes  from 
the  complexity  of  general  multibody  dynamics  systems  and  the  wide  open  design  space 
that  covers  passive,  active  and  reactive  devices  with  nonlinear,  time-dependent  and 
timing-dependent  design  variables. 
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Topology  optimization  for  optimal  structural  design  methodology  has  received  extensive 
attention  since  Bendspe  and  Kikuchi  (1988)  as  seen  by  its  wide  application  to  many 
structural  optimization  problems  (Bendspe,  1989;  1995;  Bendspe  and  Sigmund,  2003;  Ma 
et  al.,  1995b;  1995c;  Sigmund,  2001).  There  are  two  major  approaches  towards  topology 
optimization:  one  is  the  continuum  based  approach,  while  the  second  is  the  discrete 
component  based  approach.  In  the  continuum  based  approach,  the  material  is 
continuously  distributed  within  a  design  domain  by  considering  a  specific  variable 
(physical  or  artificial)  material  model  in  the  design  domain.  In  this  approach,  the 
structure  is  consequently  optimized  by  varying  the  design  variables  associated  with  the 
material  model.  In  the  discrete  component  based  approach,  for  example,  the  ground 
structure  approach  developed  by  Zhou  and  Rozvany  (1991),  a  structural  optimization 
problem  is  transformed  to  a  problem  of  seeking  the  optimal  layout  in  a  design  space  that 
considers  all  the  possible  connection  members  between  the  predefined  nodal  points  and 
the  optimization  is  achieved  by  removing  unnecessary  connection  members  and 
reinforcing  necessary  connection  members  in  the  design  space  in  improving  the  design 
objective. 

The  standard  topology  optimization  method  has  been  extended  to  a  multi-domain 
topology  optimization  (MTO)  method  (Ma  et  al.,  2006b)  to  consider  a  topology 
optimization  problem  with  multiple  domains  by  allowing  assignment  of  different 
amounts  of  the  materials,  as  well  as  of  different  materials,  to  the  different  sub-domains  of 
a  structure.  This  technique  can  be  used  to  deal  with  a  number  of  important  applications, 
such  as  structure -fixture  simultaneous  design  problems,  functionally  gradient  material 
design  problems,  and  crush  energy  management  design  problems. 
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Various  optimization  algorithms  have  been  developed  for  usage  in  topology  optimization, 
such  as  the  Optimality  Criteria  (OC)  method  by  Berke  and  Khot  (1987),  Sequential 
Linear  Programming  (SLP),  Convex  Linearization  (CONLIN)  method  by  Fleury  and 
Brainbant  (1986),  the  Method  of  Moving  Asymptotes  (MMA)  by  Svanberg  (1987), 
Diagonal  Sequential  Quadratic  Programming  (DSQP)  by  Fleury  (1987),  Modified 
Optimality  Criteria  (MOC)  method  by  Ma,  Kikuchi  and  Hagiwara  (1992)  and 
Generalized  Sequential  Approximate  Optimization  (GASO)  by  Ma  and  Kikuchi  (1995a). 
The  GASO  algorithm  extends  the  compatibility  of  previous  optimization  algorithms  by 
allowing  more  advanced  updating  rules  and  offering  more  flexibility  for  a  wide  range  of 
optimization  problems.  The  enhancement  in  the  GSAO  results  in  improved  convergence, 
higher  computational  efficiency  and  a  more  stabilized  iterative  process  for  large-scale 
optimization  problems,  including  those  dealing  with  dynamic  response.  This  method  is 
ideal  for  multi-domain  topology  optimization  problems  and  was  be  utilized  in  the  present 
effort. 

Topology  optimization  problems  usually  involve  in  a  large  number  of  design  variables; 
therefore,  an  efficient  sensitivity  analysis  method  is  critical  for  obtaining  solutions  within 
practical  time  limits.  Efficient  sensitivity  analysis  methods  have  been  developed 
previously  for  topology  optimization  related  to  static  response,  eigenvalue,  and  frequency 
response.  For  example,  Zhou  and  Rozvany  (1991),  computed  sensitivities  are  based  on 
the  static  response  of  a  linear  elastic  structural  system.  Sensitivity  calculations  for 
dynamical  systems  are,  however,  fundamentally  different  from  those  for  a  static  or  quasi¬ 
static  system.  Sensitivity  calculation  is  even  more  challenging  when  dealing  with 
multibody  dynamics  systems,  which  are  governed  by  sets  of  differential-algebraic 
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equations  (DAEs.  In  both  the  dynamic  and  multibody  dynamic  response  problems,  the 
governing  equations  are  time-dependent  and  so  are  their  sensitivities.  For  structural 
dynamic  problems,  there  are  two  widely  used  sensitivity  analysis  methods:  the  direct 
differentiation  method  and  the  adjoint  variable  method  (Hsieh  and  Arora,  1984).  To  carry 
out  sensitivity  analysis  by  the  direct  differentiation  method,  the  dynamic  equations  need 
to  be  solved  as  many  times  as  the  number  of  design  variables  (Kang,  Park  and  Arora, 
2006).  Therefore,  this  method  in  general  is  infeasible  for  topology  optimization  problems 
dealing  with  a  large  number  of  design  variables.  Cao,  et  al.  (2003)  proposed  an  adjoint 
variable  sensitivity  analysis  method  for  systems  governed  by  DAEs  of  index  up  to  two.  In 
this  approach,  a  new  set  of  DAEs  for  the  adjoint  variables  is  solved  for  obtaining  the 
sensitivities  (Alexe  and  Sandu,  2009).  For  complex  multibody  dynamics  system  models, 
the  difficulty  of  solving  the  additional  adjoint  equations  is  significant.  Recently,  for 
topology  optimization  of  a  flexible  multibody  dynamic  system,  Bruls  et  al.  (2009) 
proposed  a  sensitivity  analysis  method  based  on  the  general-a  method  (Chung  and 
Hulbert,  1993).  This  method  considers  the  dynamic  effect  of  the  multibody  dynamics 
system  based  on  the  generalized-a  method;  it  however  still  requires  solving  the  dynamic 
equations  for  each  design  variable.  Kang,  Choi  and  Park  (2001)  proposed  using 
simplified  quasi- static  load  cases  equivalent  to  the  complicated  loading  for  multibody 
dynamics  system.  However,  it  can  be  difficult  to  find  equivalent  static  loading,  and  the 
optimization  results  based  on  equivalent  static  loading  might  be  not  able  to  converge  to 
same  optimization  results  with  actual  loading  condition.  (Bruls  et  al.,  2009) 

This  paper  presents  an  extension  of  the  topology  optimization  method  for  geometrically 
nonlinear,  time-dependent  and  timing-dependent  multibody  dynamics  systems  with  the 
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consideration  of  nonlinear  response  and  a  general  multidisciplinary  system  design 
problem  with  the  various  options  from  using  passive,  active  and  reactive  structures  and 
devices.  Of  particular  emphasis  are:  1)  dealing  with  design  objectives  that  consider  time- 
and  timing-dependent,  dynamic,  large  deformation  responses;  2)  general  representative 
models  for  the  multi-disciplinary  (passive,  active  or  reactive)  components  in  a  multibody 
dynamics  simulation  system;  3)  designing  an  optimal  system  that  can  satisfy  multiple 
requirements  under  various  operating  conditions;  4)  an  efficient  sensitivity  analysis 
method  for  the  optimization  problem  of  the  occupant  restraint  system  design;  and  5)  a 
general  and  advanced  optimization  algorithm  that  can  be  used  to  solve  the  design 
problems. 

2.  Description  of  the  design  problem 

As  shown  in  Figure  1,  consider  two  general  multibody  dynamics  systems,  MDS-1  and 
MDS-2,  interconnected  by  a  set  of  N  connection  members.  Each  multibody  dynamics 
system  has  a  number  of  rigid  bodies  linked  by  joints,  bushings,  and/or  other  internal 
constraints.  As  suggested  in  Figure  1,  MDS-1  may  represent  a  human  body,  while  MDS- 
2  may  represent  a  vehicle  system.  There  are  rigid  bodies  in  MDS-1,  and  n2  rigid 
bodies  in  MDS-2.  The  set  of  connection  members  may  represent  a  possible  system  that 
restrains  the  relative  motions  between  the  two  multibody  dynamics  systems.  Each 
member  in  the  restraint  system  can  be  described  as  an  interaction  force  between  the  two 
interacting  points  at  the  two  multibody  dynamics  systems.  The  interaction  force  may 
have  non-linear  dependency  on  the  relative  movement  (displacement,  velocity,  and/or 
acceleration)  of  the  points  and  it  can  be  time-dependent  and/or  timing-dependent.  It  can 
also  be  passive,  active,  or  reactive  depending  on  the  application. 
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Figure  1  General  description  of  the  design  problem 


In  general,  the  /th  interaction  force,  which  acts  on  /nth  body  in  MDS-1  and  nth  body  in 
MDS-2,  can  be  defined  as 


./>./;(  a.  A,/.  p) 


(l) 


Here  i  denotes  the  /th  interactive  member,  A(  denotes  the  relative  distance  change 
(deformation)  between  the  two  interacting  points,  in  which  P-m)  is  the  interacting  point  of 
the  /th  interactive  member  of  the  /nth  body  in  the  MDS-1,  and  Q.n)  is  the  interacting 
point  of  the  /th  interactive  member  of  the  nth  body  in  the  MDS-2.,  A;  denotes  the  speed 
(time  directive  of  A( ),  r?  denotes  the  critical  timing  for  activating  the  /th  interactive 
member,  4°  denotes  an  initial  distance  gap  for  the  /th  interactive  member  to  become 
active,  and  p;  is  a  vector  of  other  design  parameters  for  the  /th  interactive  member.  For 
example,  a  simple  form  of  f  is  given  by: 

/  =  k  A  +  c  A.  (2) 

A  one  way  contact  with  a  gap  function  can  be  defined  as: 
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(3) 


.  _  J  0  A,  <  S' 

f'-\t:(A,-S:)  +  cA  ' 

An  active  force  function  can  be  defined  as 

/,  =  /o,  exp  (-^  (r  -  r0  )2 )  (4) 

where  ki  and  c.  are  stiffness  and  damping  coefficient  for  the  /th  interactive  member;  /0. 
and  A  are  design  parameters  for  the  /th  interactive  member. 

Since  the  /th  interactive  member  connects  the  mth  body  in  MDS-1  and  the  nth  body  in 
MDS-2,  f.  can  also  be  denoted  as  ;  A.  can  also  be  denoted  as  A'™' .  The  direction  of 

r 

the  interactive  force  of  the  /th  member  is  defined  by  e<™)  =  a  Fi  ..  ,  where  r 

. | 

denotes  the  line  of  action  between  p‘m,Q‘n> ,  Therefore,  the  /th  force  vector  acting  on  the 
MDS-1  is  f1  = ./ e. ,  and  the  force  vector  of  the  same  interaction  member  acting  on  the 
MDS-1  is  f,2  =  -ft  ,  and  we  have  f 1  +  f 2  =  0  .  Let  a  global  force  vector  F  and  global 
deformation  vector  A  be  given  as: 


II 

-  fj 

(5) 

A={A„  A2,  •• 

■  kY 

(6) 

which  represents  the  restraint  system  with  a  total  of  N  interaction  forces. 

Assume  a  global  coordinate  system  r-.o-xyz  ,  and  local  coordinate  systems 
L1,"" :  o\m) with  origin  o'”'  attached  to  the  mass  center  of  mth  body  in  MDS-1, 
4"' :  o'"' with  origin  o'"’  attached  to  the  mass  center  of  nth  body  in  MDS-2. 
Assuming  qt  =  [q(|1)r,ql2,r,...,q(,!l)r]r  is  the  generalized  coordinates  vector  of  MDS-1, 
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q2  =  [q21)r,q(12)r,...,q(2"2,r]r  is  the  generalized  coordinates  vector  of  MDS-2,  the 
governing  equation  for  MDS-1  can  be  written  as: 

M,(q,)iil-Ql+(C1)qr|>.l=F“+F1' 
c,(qrq,)  =  o 

where  the  first  equation  in  (7)  is  the  dynamic  equilibrium  equation,  and  the  second 
equation  is  the  constraint  equation  for  MDS-1.  M,  denotes  the  generalized  mass  matrix, 

(Cj  denotes  the  Jacobian  matrix  of  C, ,  denotes  vector  of  Lagrangian  multipliers.  Q, 

is  the  quadratic  velocity  term.  F,Eu  denotes  the  external  force  applied  on  MDS-1,  F,"  is  the 
generalized  force  vector  of  MDS-1  due  to  the  restraint  system  to  be  designed. 

Similarly,  the  governing  equation  for  MDS-2  can  be  written  as: 

M2(q2)q2  —  Q,  +(C,)^  X2=Ff  +  F? 

(8) 

[C2(q2,q2)  =  0 

in  which  M,  denotes  the  generalized  mass  matrix,  (C,^  denotes  the  Jacobian  matrix  of 

C2 ,  denotes  vector  of  Lagrangian  multipliers.  Q,  is  the  quadratic  velocity  term.  Ff' 
denotes  the  external  force  applied  on  MDS-2,  F2"  is  the  generalized  force  vector  of  MDS- 
2  due  to  the  restraint  system  to  be  designed. 

F®  and  F2  are  the  generalized  force  vectors  defined  in  the  generalized  coordinate  systems 
for  MDS-1  and  MDS-2.  In  general,  F®  and  F2  can  be  written  as 

F*=B'Fa„dF'=BJF  <9) 

or  equivalently, 
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=  BrF 


(10) 


F,9 

F9 

2 


where  B  =  [B,  B2J  is  called  compatibility  matrix ,  which  is  a  function  of  the 

generalized  coordinates  q,  and  q2.  B,  is  the  compatibility  matrix  for  MDS-lwhile  B,  is 
the  compatibility  matrix  for  MDS-2.  Due  to  the  nonlinear  geometry  effects,  the  B  matrix 
can  be  highly  nonlinear  with  respect  to  qj  and  q2 . 

Consider,  for  example,  a  planar  multibody  dynamics  system,  for  the  mth  body  with 
generalized  coordinates  q'"’  =  .  yo .  J  in  MDS-1,  and  the  nth  body  with 

generalized  coordinates  q'2"'  =  [x  .  yo„,  J  in  MDS-2.  Then  the  first  equation  of 

equations  (7)  and  (8)  for  the  mth  body  in  MDS-1  and  the  nth  body  in  MDS-2  can  be 
written  in  the  following  Newton-Euler  form  (Hahn,  2002): 


M[m) 

0 

0  " 

X  (M) 

oj 

LK-). 

Ar)/ 

0 

M1(m) 

0 

ya[-' 

= 

+ 

(c). 

0 

0 

j(m) 

J  1 

y/[m> 

Leir 

4i;or]-Ac,ff-[(c;-)  ( /•;; )  r 

M'm' 

Ext 

M'n) 

0 

0  " 

x ... 

o\ 

Z.M, 

in  G/‘n> 

0 

M)n) 

0 

yor 

= 

ZU*), 

+ 

(C). 

0 

0 

7(«) 

J  2 

Y?_ 

]-aA',*-[(^.4,)x  (f;j)  j' 

Ext 

(11) 


(12) 


where  M[m> ,  are  the  mass  of  the  mth  body  in  MDS-1  and  the  nth  body  in  MDS-2. 
J\""  and  J';"  are  the  moment  of  inertia  with  respect  to  mass  center  of  the  mth  body  and 
the  nth  body  respectively.  Assuming  there  are  Nm  interaction  forces  applied  on  the  mth 
body  in  MDS-1,  the  indexes  of  these  forces  elements  are  denoted  as 
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CM*' 


(m)  Am) 


/'vm)J  ,  similarly,  for  the  nth  body  in  MDS-2  we  can  define 


/<B)  =  {/'"'  ...  /wj.  Assuming  that  the  interactive  forces  apply  between  the  with  body 

in  MDS-1  and  the  nph  body,  n2th  body,  ...,  nN  th  body  in  MDS-2,  then  the  global  force 


vector  for  the  mth  body  in  MDS-1  can  be  written  as  p 


f  (™«i )  f  ( mn2  ) 

J  Am)  J  Am] 


( mnN,„  ) 


in 


which  F‘h  and  F‘h  are  generalized  forces  of  the  /  th  interactive  member  for  the  mth  body 
in  MDS-1  and  the  /.  th  interactive  member  for  the  nth  body  in  MDS-2,  expressed  in  the 

global  coordinate  system.  Note  that  x%o, . ]  and  x£,o.  J  are  the  local  position 

of  the  /.  th  attached  point  F"“  on  the  mth  body  in  MDS-1  and  the  local  position  of  the 
/..th  attached  point  Q."°  on  the  mh  body  in  MDS-2.  f)"' )  ( )  m'"/ J  and 

|j  )  m™  J  are  the  external  force  vectors  applied  on  the  respective  mth  body  in 

MDS-1  and  nth  body  in  MDS-2.  A'1  *  and  Ai  S  are  the  transformation  matrix  between 
local  coordinate  system  L'"0 ,  and  global  coordinates  system  R  . 


cos 

sin  Yi 

-  sin 

cos^;™ 

COSl/A 

sin^"’ 

-sin^,'0 

cos  y/"' 

(13) 

(14) 


The  /th  interactive  force,  which  connects  the  mth  body  in  MDS-1  and  the  nth  body  in 
MDS-2,  can  be  expressed  in  the  global  system  R  as  follows, 


1 

-“o 

_ 1 

(F,V_ 

r  r  ,  ,  (ml 

Q‘  'P'  1  Q'.'P' 

ft' 


(15) 
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Therefore,  the  global  force  vector  applied  on  the  with  body  in  MDS-1  can  be  denoted  as 

F1(",)=j f.<mni>  ...  fiian"j1^  and  calculated  as: 


<5  s 

Xb, 

i _ 

ZU‘), 

£  l-y$ x%  ]  •  AC>K  •  [( F* )  ( F* )  I' 

>m  X  y 

where 


X  1 

^  Q'"X\ 

J 

R 

rQ';i\p\l\ 

.  '■m  '>m  Jy 

“'i  'i 

Ac" 

V  '  1 

Q^P'Z]  J 

rbVmi’^o”) 

rQ{"m)P%\ 
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1  QfmT  C> 

)( 

rR„„,  1 

rs-'C; 

rn(m)  xr  / 

p\Z)°i  J 

Q(m)P(m\ 

r 

Qim)Pn  1 

/ 

> 

r 

xr  , 

/ 

r  .  , 

\(m) 

71  Q r[Z] 

Q(!trpli\ 

V  v« 

'Nm  7 

Z 

‘N„,  'Nm 

(16) 


(17) 


The  relation  between  the 


ATlth  jj(m) 

compatibility  matrix  1  in  MDS-1  and  the 


generalized 


coordinates 


Qi 


and 


are,  in  general,  highly  nonlinear. 


>«%-),] 

cos  I/S;' (r£,'o„ 

)  siny//  (/;; 

)  +sin^”)  (Xv., ) 

kv)J 

sin 

)y+cos^(^„ 

)  +  yo™  -  sin  Y\m)  ( rp-"w-> 

!y  1  \  ,  1 

)  cosy/,” )  -y 

(18) 


Substituting  equation  (18)  into  (17),  we  obtain  the  nonlinear  relation  between 
compatibility  matrix  and  the  generalized  coordinates. 

The  nonlinear  relation  between  the  deformation  of  the  hh  connecting  member  and 
the  generalized  coordinates  q|m)  =  j\v  _  v  _  and  q(,"’ =^„,1  >■  yA"’  J  is  due  to  the 


large  translation, rotation,  and  nonlinear  geometric  properties  of  dynamics  systems.  The 
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deformation  of  the  ,lh  interactive  member  attached  to  the  mth  body  in  MDS-1  and  the 
wth  body  in  MDS-2  is: 


A  ( mn )  _ 

r  „  (  ii  \ 

_ 

4»0  «  _r* 

Q-  p.’ 

r)\n  >  p\m> 

t=t0 

io 

tol 

-/° 


(19) 


Then,  the  deformation  vector  A""' for  the  mth  body  is  denoted  as 


l  ('»)  _J  a(  mnx  )  A<  mn2) 


^ mn2 


,  (mnN  ) 


r 


(20) 


The  following  relationship  is  obtained  between  the  mth  deformation  vector  A(m)  and  the 


mth  compatibility  matrix  bJ'"'  by  differentiating  equation  (20)  with  respect  to  the 


generalized  coordinates: 


d\(m) 

dqjm) 


dA(mn° 

h 

dA(mn° 

h 

dA<mn,) 

h 

dxR 

C'-A  (m) 

°\ 

°i 

ULA;(m) 

l2 

l2 

7)A(""h) 

ULA;(m) 

l2 

dxR 

O) 

°l 

dyRm) 

°l 

dy[m) 

Nm 

r)  A  (mn^m  ) 

Nm 

A  ( mnNm ) 

Nm 

dxR 

UA'0(m) 

*y% 

dy/[m) 

(21) 


3  Design  variables  in  the  optimization  problem 

The  optimization  problem  is  defined  based  on  state  equations,  general  force  elements  and 
critical  boundary  conditions.  The  design  variables  in  this  work,  a  =  [ax  a ,  aN^ , 

0<«r  <1  (i  =  l,2,...,N) ,  are  similar  to  the  relative  density  design  variables  in  power-law 
approach  or  SIMP  method,  and  are  associated  with  each  original  global  force  element  f. . 
The  design  variables  vector  a  also  could  be  defined  as  cost  functions  or  material 
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coefficients.  The  modified  global  force  element  in  the  optimization  problem  f  '  is  written 
as: 

=  (0<«,<1,  i  =  l,2,...,N)  (22) 

where  li  is  the  power  parameter 

The  global  force  vector  F  including  design  variables  will  be  rewritten  as  follows: 

F  =  [</,  </2  •••  </w]F  (23) 

4  Topology  optimization  for  multidisciplinary  structure  design 

In  general,  an  objective  function  for  multibody  dynamics  systems  can  be  written  as  a 
function  of  generalized  coordinates,  generalized  velocities  and  generalized  accelerations, 
namely,  g  =  g(q,q,q,a) .  The  topology  optimization  for  multibody  dynamics  systems  with 
multidisciplinary  structural  components  with  respect  to  dynamic  response  has  a  general 
form: 

min  g(q,q,q,a) 

a 

s.t. :  state  equations 

Z  r/a,Vi  <  h0j  {j  =  1,2,...,M)  (24) 

i=l 

0<a,  <ai  <ai  <1  (i  =  1,2,. ..,7V) 
y{  :  grouping  index  ( /'  =  0  or  1 ) 

where  M  is  the  total  number  of  constraints,  Vi  is  the  volume  or  cost  function  for  the  zth 

constraint.  The  components  in  the  restraint  system  can  be  divided  into  different  groups, 
which  may  belong  to  different  disciplines,  and  each  group  can  have  its  own  constraint, 
resulting  in  a  multi-constraint  design  problem.  Figure  2  shows  the  flow  chart  of  the 
multidisciplinary  structure  design  process. 
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Figure2  Multidisciplinary  structure  optimization  process 


4.1  Optimization  algorithm 

The  Generalized  Sequential  Approximate  Optimization  (GSAO)  developed  by  Ma  and 
Kikuchi  (1995a)  is  adopted  to  solve  this  topology  optimization  problem.  This  algorithm, 
based  on  convex  approximation,  extends  the  compatibility  of  previous  optimization 
algorithms  significantly  by  using  advanced  updating  rules  and  offering  more  appropriate 
parameters  for  the  optimization  process  algorithm.  In  specific  cases,  this  algorithm 
reduces  to  most  popular  topology  algorithms,  such  as  OC,  COLIN,  MMA,  DSPQ  and 
MOC.  The  GSAO  enhancements  result  in  improved  convergence,  higher  computational 
efficiency  and  a  more  stabile  iterative  process  for  large-scale  optimization  problems. 
GSAO  also  is  well  suited  for  multi-constraint  problems.  The  flow  chart  of  the  GASO 
optimization  process  is  shown  in  Figure  3. 

Figure3  Flow  chart  of  GSAO  optimization  process 
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Using  the  GASO  algorithm,  a  sequence  of  approximate  optimization  problems  is 
obtained: 


minimize  g*  +  ^  a.  \aj  -  cj  f 

i= 1 

i=l 

aj<ai<ai  (i-l,2,...,N) 


(25) 


By  properly  choosing  the  optimization  parameter,  the  approximate  optimization  problem 
can  always  be  made  convex.  It  is  then  solved  by  using  the  dual  method,  where  the  dual 
problem  is  given  by 

maximize  U(X  (>.),>,) 

(26) 


2,  >  0  (j  =  1,2,... ,m) 


A  typical  updating  rule  for  the  GSAO  method  is: 
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f 


Y» 


<*,  =c,.+ 


(of-c,-)  (i  =  l,2,..., AO 


v  ;=i  ; 

4.2  Sensitivity  analysis 

Combining  equations  (7),  (8),  and  (10): 


M(q)q  -  Q  +  (C)^L  =  FExt  +  BrF 

C(q,q)  =  0 


(27) 


(28) 


where, 

(C,)„  0  ■ 

0  (<y% , 

(29) 


k= 

V 

jj lExt  _ 

1 

5  ^  i 

i _ 

C  = 

c, 

>2_ 

9 

17  Lxt 

_  2 

and 

c2_ 

To  simplify  the  discussion  of  the  sensitivity  analysis  in  this  section,  it  is  assumed  that  the 
global  force  vector  F  in  equation  (28)  is  only  an  explicit  function  of  the  deformation 
vector  A  and  the  design  variables  a  ,  namely 

F  =  F(A,  a)  (30) 


q  = 


q, 

[M,  0  " 

rQii 

711 

M  = 

1 

Q= 

(c)  = 

_q2_ 

9 

1 

o 

s 

to 

1 _ 

9 

Qi_ 

v  yq 

9 

While  a  more  accurate  sensitivity  analysis  method  can  be  obtained,  we  propose  a 
simplified  but  efficient  sensitivity  analysis  method,  which  can  be  easily  implemented  into 
commercial  multibody  dynamics  codes,  such  as  MSC/ADAMS. 

The  first  equation  in  equation  (28)  can  be  rewritten  as 

F9  =M(q)q-Q  +  (C)^-F£xf  =BrF  (31) 
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Here  F''  is  the  generalized  action-reaction  force  between  the  multibody  dynamics  system 
and  the  restraint  system.  Since  the  objective  is  to  obtain  an  optimal  restraint  system,  the 
parameters  in  the  two  given  multibody  dynamics  systems  are  not  allow  to  change.  To 
apply  the  simplified  sensitivity  analysis  method,  it  is  assumed  that  Fq  =  F?(r)in  equation 
(31)  is  the  force  obtained  in  the  previous  design  stage  by  solving  equation  (28),  but  it  is  a 
given  force  when  evaluating  the  design  changes  at  the  current  stage.  This  assumption 
significantly  simplifies  the  sensitivity  analysis  process. 

Taking  the  derivative  of  equation  (31): 


0 


'dB' 

T 

'dF' 

v  da  j 

F  +  B 

v  da  j 

Similarly  from  equation  (30): 

dF  _  OF  dA  rfq  OF 
da  dA  dq  da  da 


-KB  —  +  — 
da  da 


(32) 


(33) 


and  by  application  of  the  chain  rule: 

dB  _  dB  dq 

da  dq  da 


(34) 


where  K  =  —  and  B  =  —  — . 

dA  dq 


Substituting  equations  (33)  and  (34)  into  equation  (32): 


BKB-Fr 


dB 


dqj 

which  can  be  solved  as: 


d q  =Rr  dF 


da 


da 


(35) 


dq 

da 
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BKB-F 


dB 

dq 


B 


dF 

da 


(36) 
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In  general,  assuming  objective  function  g  =  g( q,a)  is  a  function  of  generalized 
coordinates  q  and  design  variable  vector  a ,  then  we  have 


dg_=dg_dq  +  dg_=dg_ 

da  3q  da  da  3q 


^  T  T  3B 

B  KB-F 


V 


3q 


Br— +  — 

3a  3a 


Adopting  an  adjoint  vector  v ,  which  satisfies  the  following  adjoint  equation: 


BKB-Fr 

V 


3B 

v  — 

(Vf 

3qy 

V  — 

*1  ■ 

we  have 


f^L  =  vrBr  —  + 

da  da  da 


(37) 


(38) 


(39) 


For  the  special  case  where  F  =  KA  and  K  =  K(a) ,  we  will  have 


da 


Bt  —  A 
3a 


+ 


3a 


(40) 


In  general  case,  F  can  be  a  nonlinear  function  of  A  ,  but  equation  (39)  still  holds. 

4.3  Reverse  method  for  compatibility  matrix  calculation 

Generally,  the  compatibility  matrix  B  is  difficult  to  obtain,  particularly  if  the  internal 
information  of  a  multibody  dynamics  code  is  not  accessible.  There  is  a  need  to  develop  a 
more  effective  calculation  method  to  obtain  the  B  matrix  using  only  the  information 
available  during  a  normal  solution  process  without  requiring  internal  information  or 
modifying  the  multibody  dynamics  code.  In  general,  the  compatibility  matrix  B  is  the 

3A(!> 

assembly  matrix  of  the  sub-matrices  B'°  where  B(1)  =—  and  qm  is  the  generalized 

3qu; 


coordinate  vector  of  the  ith  body  in  the  multibody  system  and  A(!)  is  the  displacement 
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vector  associated  with  the  ith  body.  Assume  that  B  *' '  denotes  the  compatibility  matrix 
B(,,at  the  nth  time  step,  and  A]'1  is  the  corresponding  displacement  at  the  nth  time  step. 
Then,  using  the  first  order  Taylor  expansion  of  A^°  at  a  point  q)'1  near  toq''1 : 


dA 


(0 


Af  =  A"«  +^(q"'  -,<'>)  =  A«  -B“>  («»">-<*) 


(41) 


or 


Bli,(qli)-q^))  =  A^-A^ 


Using  the  same  process,  for  the  time  steps  n+j  (j  =  l,2,...,jn),  we  obtain 


Bin  (n(i)  -  q(,)  )  =  A(i)  -  A(,) 

Un+;  \4n+j  4o  )  ^0  ^n  +  j 


(42) 


where,  for  the  two-dimensional  system  jn  —  3  and  for  the  three-dimensional  system 


j 


n 


=  6. 


Since  A*' 1  and  qj,°  are  calculated  at  each  time  step,  by  assuming  the  compatibility  matrix 
is  constant  within  the  small  time  interval,  we  obtain 


■a(0 
4o  ’ 


q(i) 

’  ""+j„ 


-A 


(0 

n+ 1 9  *  * 


• ,  A 


(0 

0 


(43) 


Equation  (43),  can  be  solved  for  Bj,0  .  By  assembling  all  B|?°  the  global  B  matrix  is 
constructed. 


4.4  Numerical  example  for  a  two  rigid  body  dynamics  system 

A  two  rigid  body  dynamics  system  is  shown  in  Figure  4,  with  the  mass  of  body  1  =  60  kg 
and  its  mass  moment  of  inertia  =  10  kg  ■  m2 ;  and  the  mass  of  body  2  =  2,000  kg  ,  and  its 
mass  moment  of  inertia  =  1.0E6  kg  m 2 .  There  are  51  connecting  members  each  with 
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initial  linear  stiffness  =  200  N/m.  An  angular  acceleration  is  applied  to  body  2  of 


magnitude  10  rad  I  s2  with  the  rotation  center  of  02 . 


Figure  4  Two  rigid  bodies  dynamics  model 


Consider  an  objective  function  as  the  maximum  relative  translation  displacement  of  body 
1  with  respect  to  body  2,  and  the  optimization  problem  is  to  minimize  the  objective 
function.  The  objective  function  is  defined  as  follows: 


min 

a,  0=1, 2,...,N) 


maxg  = 


4 

Oi  0~) 


W 


LyR 


(44) 


where  q1 


xa  ya 


,  q2 


nr 

x0  y0  y/2\  ,  are  generalized  coordinates  body  1 


and  2,  respectively;  Au-R  is  the  transformation  matrix  from  global  coordinate  system  R 


to  local  coordinate  system  L,.  v'0\K  is  the  vector  r00  in  local  coordinate  system  at 

1  -  t—tn.  1 


the  initial  time  ,  and  W  is  a  weighting  matrix,  assumed  as  W  = 


Figure  5  Two  rigid  body  dynamics  system  optimization  result 
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Both  a  traditional  adjoint  method  and  the  proposed  sensitivity  analysis  method  were  used 
to  solve  the  example  problem.  From  Figure  5,  it  can  be  seen  that  the  adjoint  method 
converges  to  an  optimization  result  of  0.026  m2 ;  the  proposed  method  converges  to  an 
optimization  result  of  0.027  m1 .  It  is  well  known  that  computing  sensitivities  using  a 
finite  difference  method  requires  unacceptably  long  computation  times  for  a  large 
number  of  design  variables.  Using  the  adjoint  method,  it  was  necessary  to  solve  another 
set  of  differential-algebraic  equations.  The  proposed  method  calculates  the  sensitivities 
based  only  on  a  single  computation  of  the  multibody  dynamics  simulation.  Moreover 
using  the  reverse  compatibility  matrix  method  reduces  the  complexity  of  the  sensitivity 
calculation  significantly.  Therefore,  the  optimization  problem  can  be  solved  efficiently  to 
achieve  acceptable  accuracy. 

5  Application  to  vehicle  occupant  restraint  systems 

One  important  application  of  the  multidisciplinary  structure  design  methodology  is  to 
design  a  vehicle  occupant  restraint  system  to  improve  the  occupants’  safety.  The  restraint 
system  should  also  help  stabilize  the  occupant  over  rough  terrain  and  high  speed 
maneuver  conditions  and  needs  to  be  user  friendly,  such  as  easy  to  put  on  and  take  off. 
The  restraint  system  involves  a  possible  usage  of  passive,  active,  and  other  reactive 
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devices,  and  these  multidisciplinary  safety  elements  such  as  belts,  airbags  or  retractors 
have  to  be  activated  in  a  specific  sequence  and  timing  to  protect  the  occupant  in  extreme 
conditions.  Minimizing  the  system  weight,  cost  and  complexity  are  also  considered  in 
design  process.  Therefore,  it  is  necessary  to  develop  a  general  and  systematic  design 
approach  and  optimization  tool,  which  can  enlarge  the  design  space  and  obtain  optimal 
layout  design  for  best  performance/weight  and  performance/cost  ratios.  Traditional 
design  solution  based  on  engineers’  intuition  may  not  provide  the  best  combination  of 
functionality. 

Virtual  prototyping  multibody  dynamics  models  are  developed  for  computational 


simulation  in  a  commercial  code.  The  detailed  specification  of  a  virtual  24-years  old  male 


occupant  multibody  dynamics  model  (Figure  6)  is  listed  in  Table  1. 


Table  1  Specifications  of  occupant  model 
Weight 

77  Kg 

Height 

1.778  m 

CGX  (  +:  rearward  from  the  front  axial) 

1.848  m 

CGY  (+:  rightward  from  midplane) 

0.041  m 

CGZ  (+:  upward  from  the  ground) 

1.758  m 

Part  number 

58 

There  are  three  connecting  bushings  created  for  integrating  the  occupant  and  vehicle 


model  together  (Figure  7).  Two  bushings  connect  the  occupant’s  hands  with  the  vehicle, 


and  one  bushing  connects  the  occupant’s  lower  torso  with  the  seat  on  the  vehicle  to 


simulate  the  occupant’s  sitting  posture.  The  detailed  specifications  of  the  integrated 


model  are  given  in  Table  2.  The  joint  stiffness  properties  of  the  occupant  are  based  on  the 


data  measured  from  a  Hybrid  III  dummy  finite  element  model  in  a  software  library  and 


biomechanical  publications  (Dhaher  et  al.,  2005,  Dinant  and  Kistemaker,  2007;  Granata 
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et  al.,  2004;  Gunther  and  Blickhan,  2002;  Leger  and  Milner,  2000;  LSTC,  2007; 
Magnusson,  1988,  Van  der  Spek  et  al.,  2003;  Xu,  1999). 

Figure  7  Integrated  occupant  and  vehicle  model  Table  2  Specifications  of  integrated  model 


Weight 

2898  Kg 

CGX  (  +:  rearward  from  the  front  axial) 

1.70  m 

CGY  (+:  rightward  from  midplane) 

0.0006  m 

CGZ  (+:  upward  from  the  ground) 

0.80  m 

Part  number 

128 

Three  virtual  proving  grounds  were  employed  for  this  study:  severe  braking,  rollover  and 
rough  terrains.  For  the  severe  braking  case,  the  initial  vehicle  longitudinal  velocity  is  17 
m/s,  and  the  vehicle  deceleration  is  9.8  m/s  .  For  the  rollover  case,  the  initial  vehicle 
longitudinal  velocity  is  17  m/s,  and  the  steering  wheel  rotates  720°  in  1  second  for  the 
vehicle  system.  For  the  rough  terrain  case,  the  initial  vehicle  longitudinal  velocity  is  17 
m/s,  and  the  road  profile  is  a  sinusoid  function  with  magnitude  of  0.05  m  and  wave  length 
of  8  m. 

Critical  conditions  for  the  design  problem  were  identified  and  design  uncertainties  were 
eliminated,  including  the  effect  of  the  gunner’s  awareness  in  terms  of  hand  grasping:  i) 
gunner  intentionally  grasps  the  handle  in  a  maneuver;  ii)and  gunner  does  not  grip  handle 
in  a  maneuver;  the  effect  of  hand  gripping  strength  with  a  stronger  gunner  and  a  weaker 
occupant;  the  effect  of  joint  stiffness  where  the  gunner  intentionally  holds  the  position  or 
the  gunner  is  in  the  relaxed  condition;  the  effect  of  terrain  roughness  with  rough  terrain 
and  flatter  terrain;  the  effect  of  gunner  postures  considering  seated  and  standing  postures 
with  different  orientations.  Detailed  results  can  be  found  in  Dong  et  al.  (2009)  and  Ma  et 
al.  (2010).  As  an  example  of  these  studies  to  identify  critical  conditions,  consider  of  the 
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gunner’s  center  of  gravity  (CG)  height  with  respect  to  the  vehicle’s  roof,  as  shown  in 


Figure  8  and  Figure  9.  The  gunner  is  ejected  during  severe  braking  if  the  gunner’s  hands 
are  not  grasping  anything  on  the  vehicle,  but  remains  in  the  crew  compartment  in  hands 
grasping  case.  It  is  concluded  that  the  condition  of  hands  free  grasping  is  more  critical  in 
the  restraint  system  design. 

Figure  8  Occupant  CG  relative  height  response  Figure  9  Occupant  response  in  brake  condition 
with  different  grasping  condition  in  brake  at  3s  (a)  hand  free  (b)  hand  grasping 
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The  initial  design  space  was  set  up  with  evenly  distributed  connecting  members  with 
linear  stiffness  in  all  possible  connections.  In  order  to  discretize  the  design  space,  as 
shown  in  Figure  10,  5  vertical  layers  of  connection  nodes  were  placed  on  the  vehicle,  22 
predetermined  connection  nodes  on  the  occupant,  resulting  in  580  connection  members 
between  the  gunner  and  the  vehicle.  The  function-oriented  multidisciplinary  structure 
optimization  was  employed  to  optimize  the  geometrically  nonlinear,  time-dependent 
structural/multibody  dynamics  system  based  on  the  connectivity  of  interaction  points  on 
occupant  and  vehicle,  optimal  interaction  members  and  optimal  physical  properties  of  the 
interaction  members  for  occupant  at  vehicle.  The  optimal  structure  layout  was  obtained 
by  removing  unnecessary  connecting  members  and  reinforcing  necessary  connecting 
members  via  the  optimization  algorithm.  Critical  for  the  optimization  was  the  use  of  the 
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proposed  sensitivity  calculation  to  efficiently  address  the  nonlinear  geometry  effects  and 
large  motion  in  dynamic  response  optimization. 

Figure  10  Initial  structural  universe  with  connecting  points  on  (a)  vehicle  and  (b)  occupant 


Every  connecting  member  was  associated  with  a  material  coefficient,  0  <  OCi  <1  for  the 
ith  member,  and  it  was  assumed  that  stiffness  is  proportional  with  this  material 
coefficient,  so  the  stiffness  assigned  for  the  ith  connecting  member  is  cxk{) .  Consequently, 
a.  =  0  means  this  member  should  be  removed  in  the  layout,  and  a,  =  1  means  this 
member  should  remain. 

The  flow  chart  in  Figure  11  shows  the  optimization  process  applied  for  an  occupant 
restraint  system  optimization  design,  which  was  implemented  by  coupling  the 
commercial  codes. 

Figure  11  Flow  chart  of  multidisciplinary  structure  optimization  process 
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In  order  to  avoid  an  physically  infeasible  and  expensive  restraint  system,  the  design  space 
was  reduced  to  only  keeping  180  connecting  members  between  the  gunner’s  upper  torso, 
central  torso,  lower  torso  and  vehicle  body  for  the  preliminary  study.  The  maneuver 
condition  for  the  vehicle  is  a  step  steer  condition,  in  which  the  steering  wheel  rotates  360° 
in  0.5  seconds  with  initial  longitudinal  velocity  of  17m/s.  The  objective  function  is 
defined  as  the  maximum  relative  translation  displacement  between  the  gunner’s  center  of 
gravity  and  the  vehicle  in  the  time  duration  [?„,?,] ,  and  design  objective  is  to  minimize  the 
maximum  relative  translation  displacement  as  defined  in  equation  (45) 

C  r  -1 T  r  -d 


min  <:  max  g  = 

<Xt (i=l,2,..„ N)  |  [f0l/f 


A',*(q0-qv)-r1 


OnOv 


w 


A'-'V-qO-r, 


(45) 


where  q° , ql  are  the  respective  generalized  coordinates  for  the  gunner’s  center  of  gravity 
and  vehicle  body,  and  W  =  ^ 

Figure  12  shows  the  design  objective.  From  the  design  objective  iteration  results  in 
Figure  12  (a),  it  is  concluded  that  the  proposed  function-oriented  design  method  based  on 
topology  optimization  can  solved  the  problem  appropriately  and  also  reduce  the 
maximum  occupant  relative  translation  displacement  with  respect  to  the  vehicle  with 
fewer  active  connecting  members.  That  is,  the  active  remaining  members  can  restrain  the 
occupant  at  the  initial  position  more  effectively  than  the  initial  evenly  distributed 
members.  In  Figure  12  (b),  the  black  color  denotes  a  higher  value  of  design  variable,  i.e., 
the  member  should  remain  in  the  final  layout;  while  the  grey  color  means  a  medium 
value  of  design  variable,  or  the  members  need  more  investigation  and  the  white  color 
means  a  lower  value  of  design  variable,  or  the  members  can  be  removed  in  the  final 
layout. 
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Figure  12.  Optimization  iteration  results  (a)  and  interactive  members’  final  stiffness 
distribution  (b) 


Figure  13  shows  the  gunner’s  dynamic  response  under  the  step  steer  condition.  The  upper 
row  is  the  dynamic  response  before  optimization,  in  which  all  the  interactive  members’ 
stiffness  was  distributed  evenly;  the  lower  row  is  after  optimization,  with  stiffness 
distributed  as  shown  in  Figure  12  (b).  From  the  rightmost  frames,  it  is  obvious  that  the 
optimized  interactive  members  layout  can  constrain  the  occupant  much  more  effectively 


in  dynamic  loading  condition  with  same  total  stiffness  amount. 


Figure  13.  Occupant  dynamic  response  in  (a)0.25s  (b)0.5s  (c)0.75s  and  (d)  Is 


(a)  0.25s  (b)U.Ss  (cjlL/hs  (d)  is 

Further  interpretation  of  the  optimum  layout  depends  on  the  mechanical  properties  of  the 


remaining  connecting  members  and  the  engineer’s  intuition.  If  deformation  of  remaining 


members  is  further  investigated,  the  compression  members,  which  are  shown  as  (-),  can 
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be  realized  as  airbag  devices  and  the  tension  members,  which  are  shown  as  (+),  could  be 


realized  as  belt  devices  in  further  components  design  based  on  Figure  14. 
Figure  14.  Optimum  connecting  members  layout 


The  ordinary  linear  stiffness  members  in  the  multibody  dynamics  system  can  also  be 
substituted  with  connecting  members  with  nonlinear  stiffness,  such  as  shown  in  Figure 
15.  Using  this  nonlinear  stiffness  response,  the  optimized  result  under  the  step  steer 
condition  is  shown  in  Figure  16,  and  it  is  concluded  that  the  proposed  function-oriented 
design  method  based  on  topology  optimization  can  also  be  applied  to  the  system  design 
with  nonlinear  interactive  members. 

Figure  15.  Interactive  members’  nonlinear  Figure  16.  Nonlinear  members’ 
stiffness  optimization  results 
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6  Critical  parameters  identification  for  multidisciplinary  components 
design 

Once  obtaining  an  optimum  layout  using  the  proposed  general  multidisciplinary  structure 
function-oriented  design  method,  it  is  necessary  to  identify  the  type  of  interactive 
members,  such  as  passive,  active,  reactive,  and  to  identify  timing  parameters,  length  of 
actuation  period  or  other  effective  design  variables,  such  as  critical  design  parameters  of 
mechanical  properties  for  different  components.  To  do  so,  various  nonlinear  general  force 
(G-force)  elements  need  to  be  developed  to  represent  multidisciplinary  components,  and 
then  be  incorporated  into  the  design  problem. 

As  an  example,  consider  a  belt  retractor  design,  in  which  the  number  of  retractors,  single 
point  or  multiple  points;  location  of  the  connecting  points,  both  on  the  occupant  and  on 
the  vehicle  are  obtained  in  the  final  optimum  layout  design  using  the  proposed  method. 
Critical  design  parameters  for  properties  of  the  retractor  need  to  be  identified  in  next  step. 
A  series  of  five  bench-top  retractor  tests  were  conducted  by  Newberry  et  al.  (2006)  on  a 
typical  pyrotechnic  retractor.  Based  on  the  experimental  data  of  Newberry,  a  G-force 
element  for  the  retractor  is  given  as 

F(t)  =  F0e~ao('~ro)1  (46) 

The  critical  design  parameters  for  the  retractor  include  peak  timing  t0 ,  pulse  width  <r0 
and  peak  value  F0.  As  a  first  comparison,  2000N  is  chosen  as  base  line  of  peak  force 
value.  3.0s/3.1s/3.2s  are  selected  as  different  peak  timing  to  investigate  the  peak  timing 
effect.  The  maneuver  condition  is  the  rollover  case.  As  shown  in  Figure  17,  it  can  be  seen 
that  later  peak  timing  causes  a  higher  possibility  of  occupant  ejection.  Peak  timing  is 
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critical  for  the  design  because  earlier  peak  timing  can  be  difficult  to  determine  by 


sensors  assessing  whether  or  not  rollover  may  happen,  and  later  peak  timing  may  not  pull 
the  occupant  into  the  crew  compartment. 

Figure  17  Comparison  for  retractor  peak  Figure  18  Comparison  for  retractor  peak 
timing  force  value 


Time  (s) 


Time  (s) 


As  a  second  study,  3.0s  is  chosen  as  the  base  line  of  retractor  peak  timing,  and 
1000N/1500N/2000N  are  selected  as  different  peak  force  values  to  investigate  the  peak 
force  effect.  From  the  results  in  Figure  18,  it  is  seen  that  smaller  peak  force  values  cause 
a  higher  possibility  of  occupant  ejection.  The  peak  force  is  critical  for  critical  for  the 
design  because  smaller  peak  force  value  could  not  pull  the  occupant  into  the  crew 
compartment,  however,  a  larger  peak  force  value  has  more  possibility  to  cause  injury  to 
the  occupant. 

It  is  concluded  that  peak  timing  and  peak  force  value  as  parameters  of  retractor  property 
are  critical  for  the  component  design.  A  representative  general  force  element  for  the 


retractor  should  include  these  two  parameters  as  design  variables. 
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7.  Conclusion 


Fundamental  multidisciplinary  structure  design  technology  is  proposed  for  a  multibody 
dynamics  systems  design  problem  which  may  have  various  options  associated  with  using 
passive,  active,  and  reactive  devices  or  materials. 

The  proposed  optimization  design  method  can  deal  with  objective  functions  that  are 
related  to  dynamic  responses  of  multibody  dynamics  systems  rather  than  static  response, 
and  that  satisfy  multiple  requirements,  such  as  those  for  designing  a  vehicle  occupant 
restraint  system,  under  various  operating  conditions  and  performance  requirements.  The 
proposed  advanced  topology  optimization  technique  uses  an  efficient  sensitivity  analysis 
technique  necessary  for  practical  multidisciplinary  multi-constraint  problems.  A 
representative  model  for  multidisciplinary  components,  including  possibly  passive,  active 
and  reactive  devices  was  developed  for  identifying  an  optimal  layout  from  a  wide  open 
design  space. 

Disclaimer 

Reference  herein  to  any  specific  commercial  company,  product,  process  or  service  by 
trade  name,  trademark,  manufacturer,  or  otherwise,  does  not  necessarily  constitute  or 
imply  its  endorsement,  recommendation,  or  favoring  by  the  United  States  Government  or 
the  Department  of  the  Army  (DoA).  The  opinions  of  the  authors  expressed  herein  do  not 
necessarily  state  or  reflect  those  of  the  United  States  Government  or  the  DoA,  and  shall 
not  be  used  for  advertising  or  product  endorsement  purposes. 
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